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ABSTRACT 

The jet of the radio galaxy M87 is misaligned, resulting in a Doppler factor 6 ~ 1 for 
emission of plasma moving parallel to the jet. This makes the observed fast TeV flares on 
timescales of f v ~ 5R g /c harder to understand as emission from the jet. In previous work, we 
have proposed a jets-in-a-jet model for the ultra-fast TeV flares with f v <c R g /c seen in Mrk 
501 and PKS 2155-304. Here, we show that about half of the minijets beam their emission 
outside the jet cone. Minijets emitting off the jet axis result in rapidly evolving TeV (and 
maybe lower energy) flares that can be observed in nearby radio galaxies. The TeV flaring 
from M87 fits well into this picture, if M87 is a misaligned blazar. 

Key words: galaxies: active - galaxies: individual: M87 - galaxies: jets - radiation mecha- 
nisms: non-thermal - gamma rays: theory 



1 INTRODUCTION 

The jet originating from the nucleus of M87 has been well stud- 
ied at all wavelengths thanks to its proximity. Its radio images and 
modeling of its interaction with its environment suggest that the jet 
is misaligned by 9 ~ 30° with respect to our line of sight (Biretta, 
Zhou & Owen 1995; Bicknell & Begelman 1996). 

Recently several groups (Aharonian et al. 2006; Albert et al. 
2008; Acciari et al. 2008) have reported TeV emission from M87 
that shows rapid variability on timescales of ~ 1 day. The fast flar- 
ing implies emitting regions of typical length scale of / cm < 5 x R g S, 
where R g = GM/c 2 and M = 3 x 10 9 M G is the mass of the black 
hole measured by gas dynamics in the vicinity of the black hole 
(Macchetto et al. 1997 j_J and S is the Doppler factor of the emit- 
ting material. Since the jet from M87 is misaligned by ~ 30°, any 
Doppler factor associated with the bulk motion of the jet is S ~ 1, 
leading to / em < 5R g . 

Aharonian et al. (2006) proposed that the TeV emitting region 
is not associated with the jet but with the magnetosphere of the 
black hole (Levinson 2000). Alternatively, the TeV emission may 
originate in the sub-pc jet while the jet is collimating (Lenain et al. 
2008) or in (at least) two zone models such as a jet with a fast spine 
and a slower layer (Tavecchio & Ghisellini 2008) or a decelerating 
jet (Georganopoulos, Peiiman & Kazanas 2005). The knot HST-1 
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1 Gebhardt & Thomas (2009) claim that the black hole mass is a factor of 2 
larger than previously thought. While in this paper, we keep M = 3x 10 9 M o 
as the reference value for the black hole mass in M87, a larger black hole 
corresponds to more extreme variability of the TeV emission that strength- 
ens our main arguments. 



has also been discussed as a possible source of the TeV emission 
(Cheung, Harris & Stawarz 2007). 

The latest (February 2008) observed TeV flaring of M87 was 
caught simultaneously in radio VLBI and X-rays (Acciari et al. 
2009). At the period of the TeV activity, the radio started to show a 
gradual rise that continued over the succeeding months. The VLBI 
map showed new radio blobs starting to move outward from within 
~ 100/?„ from the black hole, suggesting that the gamma-ray flare 
comes from very close to the black hole. While the X-rays from the 
nucleus were climbing during the flare, HST- 1 was in a low X-ray 
state, making it unlikely that HST-1 is the source of the flares. 

Fast flaring on timescales f v ~ R g /c observed in blazars is 
attributed to large Doppler factors S » 1 of emitting plasma asso- 
ciated with relativistic jets pointing directly at us. The very fast TeV 
variability of the blazars Mrk 501 and PKS 2155-304 (Aharonian et 
al. 2007; Albert et al. 2007) corresponds to f v ~ 0AR g /c and poses 
strong constraints to any model. For the TeV radiation to avoid ab- 
sorption at the source, the emitting material must move with a bulk 
Lorentz factor r cm > 50 (Begelman, Fabian & Rees 2008) while 
resolved patterns on sub-pc scales are compatible with moderate 
Lorentz factors Tj < 10 for the jet (Piner & Edwards 2004; Giro- 
letti et al. 2004). In a previous work, we proposed that very rapid 
TeV variablity observed in blazars comes from blobs of energetic 
particles that are moving with relativistic speed relative to the mean 
flow of the jet (Giannios, Uzdensky & Begelman 2009). We asso- 
ciate these minijets with the material outflowing from reconnection 
events in a Poynting-flux dominated jet. lets with typical Tj < 10 
(as indicated by proper motions) would thus contain minijets that 
move with r cm ~ 100 » Tj relative to the observer, as needed for 
TeV radiation to escape the source and power flares with t r <k R g /c. 

A natural consequence of the jets-in-a-jet model is that, while 
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the emission of about half of the minijets is beamed at an angle 9 < 
1 /Tj with respect to the jet axis, the rest emit outside the emission 
cone of the jet. In this paper, we explore the emission from these 
"misaligned" minijets. We show that the probability that a minijet 
is observed to emit at an angle 9 depends rather weakly on 9. We 
calculate the energetics and timescales of emission from minijets 
pointing off the jet axis. Off-axis minijets are bright enough to be 
observed in nearby radio galaxies. Finally, we apply the jets-in-a-jet 
model to the fast TeV flaring observed in M87. 



2 JETS IN A JET 

Our main idea of "jets-in-a-jet" is that dissipation of magnetic en- 
ergy in the jet (e.g., the result of jet instabilities) leads to minijets 
of energetic particles that move relativistically within the jet. Emis- 
sion from the minijet results in the observed short-timescale flaring. 
Depending on the direction of motion of the minijet (in the jet rest 
frame), a minijet emits within the cone of the jet or outside it. In 
jets moving toward our line of sight (blazars), emission from mini- 
jets results in powerful, rapidly evolving flares while in misaligned 
jets (e.g., in radio galaxies), minijets lead to weaker but still rather 
rapid flares. 

We believe that jets-in-a-jet may be expected in a Poynting- 
flux dominated flow. They can form in reconnecting regions of 
strongly magnetized plasma. The magnetic reconnection may result 
from the nonlinear evolution of kink instabilities in the jet (Eich- 
ler 1993; Begelman 1998; Appl, Lery & Baty 2000; Giannios & 
Spruit 2006; Moll, Spruit & Obergaulinger 2008; Moll 2009; but 
See McKinney & Blandford 2009) or reversals of the magnetic field 
polarity in the inner disk/black hole magnetosphere (Giannios et al. 
2009). Interestingly, recent evidence for such reversals come from 
the observed "flip" in the gradient of the Faraday rotation across 
the jet of B 1803+784 between the 2000 and 2002 observations of 
the source (Mahmud, Gabuzda & Bezrukovs 2009). A similar sce- 
nario of relativistic motions within the jet has also been explored in 
the context of gamma- ray bursts (GRBs). Minijets (caused by rel- 
ativistic turbulence or magnetic reconnection) may be responsible 
for the observed variability of the GRB emission (Blandford 2002; 
Lyutikov 2006a,b; Narayan & Kumar 2008; Lazar, Nakar & Piran 
2009). 

Consider a jet that moves radially with bulk Tj, in which a 
minijet develops with a Lorentz factor r co at an angle ff with re- 
spect to the radial direction (both measured in the jet rest frame). 
All primed/tilded quantities are measured in the rest frame of the 
jet/minijet respectively. In the lab frame the minijet moves with 

Fem — rjr co (i + VjV co COS0') (1) 

and at an angle 

V™ sin ff 

tan 9 = (2) 

rj(v co cos 9' + vj) 

with respect to the radial direction. For ff ~ n/2 and v co ,Vj ^ 1, 
expressions (1) and (2) give r cra ~ FjF co and 9 ~ l/Ts. Minijets 
moving at 6 > 1 /Tj have more moderate bulk Lorentz factor in the 
lab frame but, in many cases, larger than that of the bulk of the jet 
Tj (see next section). 

From eq. (2), one can immediately see that for ff < n/2 the 
emission takes place at 9 < \/T- p i.e., within the emission cone 
of the jet and vice versa. In the magnetic reconnection picture, ev- 
ery dissipation event results in a pair of mini jets that leave the 
reconnection region at opposite directions 9' and n — ff. Since the 



distribution of 6' of the minijets is symmetric around ff = n/2, for 
every minijet that is beamed within the jet 1 /Tj emission cone there 
is another minijet that points outside. 

2.1 Distribution of emitting angles of minijets with respect to 
the jet axis 

In this section, we explore how the angular dependence of the emis- 
sion from the minijets is connected to their assumed distribution of 
directions in the jet rest frame. We find that the probability that the 
observer is located within the emission cone of a minijet is rather 
insensitive to the inclination of the jet with respect to the observer. 

We consider a relativistic conical jet with opening angle 9 S ~ 
1 /T s . Emission associated with the bulk motion of the jet is beamed 
at the observer when the latter is located at an angle 9 < 1 /Tj with 
respect to the jet axis. An "off-axis observer" is one for which 6 > 
\/Tj. We measure how "off-axis" the observer is located by the 
parameter a defined as: 

a 

e=-- O) 

1 i 

If an observer is located at an angle 6 with respect to the jet 
axis, how likely is it that the emission from a minijet takes place 
within the observer's cone? As discussed above, in the lab frame, 
the majority of minijets point at small angle 9 or corresponding 
a< a few. On the other hand, minijets pointing more off-axis have 
lower bulk Lorentz factors and therefore, wider beaming angles as 
we will see. The net result is that there is almost uniform probability 
that a minijet is observed at any angle 9. 

To quantify this statement, we consider several different cases 
for the distribution P(Q.')dQ.' of the angle ff of the minijets in the 
jet rest frame. We consider distributions that are symmetric around 
ff = n/2, as expected from the formation of pairs of counter- 
streaming minijets that result from magnetic reconnection. We con- 
sider 3 cases: P(Q!) oc sin" 9', for n = 0, 1,2. The n = case cor- 
responds to the isotropic distribution. For n = 1,2 the distributions 
increasingly peak around 9' = n/2. We think it is very plausible 
that the reconnecting current sheets are preferentially oriented per- 
pendicular to the jet direction (e.g., a sequence of toroidal fields 
with reversing radial currents) leading to the ejection of minijets 
preferably at 9' ~ n/2. 

We assemble a large number of minijets with & drawn from 
the above distributions and calculate the distribution P(a)d(a) of 
emission angle a using the expression (2). For the results shown in 
Fig. 1, we have set Tj = r co = 10. Very similar angular distribu- 
tions are found for other values of Tj, r co » 1. The distribution of 
direction of the minijets peaks at a ~ 1 (i.e., at angle 9 ~ \/T } ). 
For increasing n, the 0.5 < a < 2 angles become more likely. This is 
expected since the number of minijets moving at ff ~ n/2 (which 
results in a ~ 1) increases with n. Furthermore, by definition of the 
symmetry of the distribution around ff = n/2, about half of the jets 
point at a > 1, independently of n. 

To calculate the probability distribution that a blob is observed 
at an angle 9 = a/Tj one has to take into account that the emission 
of a blob is beamed into a narrow cone of fi cm oc l/r^ m . In Fig. 2, 
we show the distribution P(cr)dSl/Tl m which we refer as the "ob- 
served angular distribution" P obs (a)dfl of the minijets. 

From Fig. 2, it is clear that for n = the observed angular dis- 
tribution of the minijets is almost isotropic. For n = 1, 2 it slightly 
favors 0.5 < a < 3. The reason for P obs being much flatter than P is 
that, although most minijets (per unit of solid angle) point within 
the jet cone, their emission is beamed within a very narrow cone. 
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Figure 1. Distribution P(a)da of direction of minijets in the lab frame 
(where a = 0T,) for different distributions P(n')dfl' of direction of the 
minijets (in the jet rest frame). 



Figure 3. Fractional accuracy of approximate expression |4) of the bulk 
Lorentz factor of the minijet compared to the exact expressions (1) and (2) 
as function of inclination a. The different curves correspond to different 
values of Tj and r co . For a>2, the approximate expression gives the correct 
answer within better than 20%. 
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Figure 2. Probability (not normalized) P obs (a)dCl oc P(a)d(Di)/r% m that a 
minijet is observed at angle a for different distributions P(£l' )dfi' of direc- 
tion of the minijets (in the jet rest frame). 



Minijets pointing to larger angles are more broadly beamed. As a 
result, the probability that a blob is observed at an angle depends 
weakly on 9. 

Out of N minijets forming in the jet, approximately N/Tl a will 
be beamed at an observer located at any angle with respect to the 
jet axis. The observed duration and intensity of these flares will, 
however, depend on the observer's angle with respect to the jet axis, 
with on-axis observers seeing far more powerful flares. We quantify 
this point in the next section. 



2.2 Dependence of observed properties of the minijets on the 
inclination to the jet axis 

As we have shown in the previous section, the probability that 
a minijet emits within the observer's line of sight depends only 
weakly on the inclination of the observer with respect to the jet axis. 
However, the detectability of emission of a minijet depends on the 
luminosity of the source. The bolometric luminosity of a minijet is 
related to the total radiated energy of the minijet E, the solid angle 



Q em over which the emission takes place, and the observed duration 
6t obs of the emission: L = E/Cl cm St obi . In the following we discuss 
how E, fi em and 6t obs depend on the inclination of the minijet. 

The expressions (1) and (2) give the velocity of the minijet 
in the lab frame once the comoving angle 8' is specified. Here, 
we derive much simpler, approximate expressions for the veloc- 
ity and direction of motion of the minijet that points at 6 > l/Tj 
(off the jet axis). In the following we express the various quantities 
as functions of the inclination a. In the limit where Tj, r co » 1 and 
1 ~ oi 5 Iji r c o, eq. (1) for the Lorentz factor of the blob becomes 



2r T 



and eq. (2) 
2 



(4) 



(5) 



T s (7i -6') n- 0' 

In Fig. 3, we show the fractional error of expression © with 
respect to the exact one [eq. (1)] as function of inclination a and 
for different values of Tj, r co . For T:, r co > 5, the approximate ex- 
pression is rather accurate provided that a > 2. In the analytical 
estimates that follow we use the approximate expression 10. 

From eq. (4), we find that a minijet moves relativistically to- 
wards the observer for a wide range of inclinations a < (2Fjr co ) 1/2 . 
This leads to the possibility of emitters with high Doppler factor 
5 » 1 even in cases where the jet is misaligned by large angles 
6»>30°. 

We assume that all the minijets are identical in a frame co- 
moving with the jet and differ only in their orientation 8'. Sup- 
pose that the energy of the minijets is E (in the minijet rest frame). 
When 6 < l/Tj (i.e., a < 1), a blob has lab-frame energy E on = 
T em E ~ Fjr c0 £ while for > l/Tj the energy of the minijet is 
Egg ~ Tr^T^E/a 2 ~ 2E OI[ /a 2 . Since the emission of the minijet 
takes place within a solid angle f2 cm ~ l/r^ m , the emission from 
minijets that point off the jet axis is less beamed than that of on- 
axis ones: Q. off ~ a 4 fl on /4. 

The observed timescale of emission from a minijet depends 
on the size of the blob, its bulk Lorentz factor, and the radia- 
tive mechanism or the duration of the reconnection event. Identi- 
cal blobs (i.e., of the same size, density, magnetic field strength) 
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that emit through the synchrotron and synchrotron-self-Compton 
(SSC) mechanisms have the same emission timescale St in the rest 
frame of the minijet. The observed duration of the blob emission 
is dt obs = St /T,. m . Thus, the duration of the flare for an off axis ob- 
server ijl^foff ~ a 2 St OI J2. This expression shows that the timescale 
of emission becomes longer when the jet is viewed from larger in- 
clination. On the other hand, the dependence on a is rather mod- 
erate. If minijets result in the very rapidly evolving flares of the 
blazars Mrk 501 and PKS 2155-304 with dt on ~ 0AR g /c (Giannios 
et al. 2009), misaligned jets with a ~ a few can also lead to rather 
short variability timescales St e ~ R g /c. 

Combining the above estimates, if the minijets emit through 
synchrotron and SSC, the bolometric luminosity during the flare of 
an off-axis blob is 



Although the a- 8 dependence is very steep, making flares from mis- 
aligned jets rather weak, it is still possible to observe them from 
nearby sources. 

The gamma-ray emission from the minijets may not be dom- 
inated by SSC, but by inverse Compton scattering of an external 
source of seed photons, provided that a powerful enough source 
exists. The accretion disk and the broad line region have been fre- 
quently considered in this respect (Dermer, Schlickeiser & Mas- 
tichiadis (1992); Sikora, Begelman & Rees (1994)). In the jets- 
in-a-jet model discussed here, soft (synchrotron) emission emitted 
at a different location in the jet can also be considered as an ex- 
ternal photon field from the point of view of the minijet. Such a 
photon field, approximately isotropic in the rest frame of the jet, 
may be connected to the process responsible for the quiescent jet 
emission or may come from the averaged synchrotron emission of 
all the minijets forming during the flaring activity of the jet. One 
can show that the expression (6) for the dependence of the blob's 
luminosity on its inclination to the jet axis holds for the external 
inverse Compton mechanism, provided that the photon field is ap- 
proximately isotropic in the rest frame of the jet. On the other hand, 
if there is a significant source of soft radiation coming from the base 
of the jet (e.g., the accretion disk), the scaling (6) fails. 

As we shall see in the next section, the observed soft emission 
coming from the nucleus of M87 (associated with the accretion disk 
and/or the inner jet) provides plenty of seed photons to be inverse 
Compton scattered in the minijet. External inverse Compton ap- 
pears to be an important mechanism for the gamma-ray emission 
from this source. 



3 APPLICATION TO THE FLARING TEV EMISSION 
FROM M87 

The nucleus of M87 contains a M = 3 x 10 9 M Q black hole (corre- 
sponding to light-crossing time of R g /c ~ 1.5 x 10 4 sec) measured 
from gas kinematics on scales of 10's of parsecs (Macchetto et al. 
1997). The M87 jet, believed to originate from the black hole, has 
been extensively observed from radio to TeV and from sub-pc to 
kpc scales. 

Modeling of the kpc-scale interaction of the jet with the ex- 
ternal medium constrains the properties of the jet. Though rather 

2 If the observed timescale of emission is determined by the duration of 
the reconnection event, the duration of the flare for an off-axis observer can 
similarly shown to be St a ff ~ a 2 St an /2. 



model dependent, the jet inclination is estimated to be 9 ~ 30° with 
respect to the line of sight, the bulk Lorentz factor of the jet Tj ~ 5 
and the kinetic luminosity Lj itme ~ 10 43 erg/sec (Bicknell & Begel- 
man 1996; Reynolds et al. 1996, Owen et al. 2000; Stawarz et al. 
2006; Bromberg & Levinson 2009). 

The nucleus of M87 and the knot HST-1 (about 60 pc away 
from the nucleus) are active with luminosity of order of ~ 10 41 
erg/sec each over a wide range of wavelengths (from IR to gamma- 
rays; see, e.g., Sparks, Biretta & Macchetto 1996; Perlman et al. 
2001). Both the nucleus and the knot HST-1 have shown variability 
in X-rays and optical on timescales ranging from months to years 
(see, e.g., Harris et al. 2009) indicating rather compact emitting 
regions. During 2005, M87 was luminous in the TeV band with 
isotropic Lx c v ~ 3 x 10 40 erg/sec and exhibited flares evolving on 
~ 2-day timescales (Aharonian et al. 2006). The poorer angular 
resolution in TeV energies allows for both the nucleus and HST-1 
as possible sources of the high-energy emission. During the 2005 
TeV flaring, HST-1 was at its peak of its optical and X-ray activity 
(exceeding the nuclear emission), making the knot a promising can- 
didate for the origin of the flares (Cheung et al. 2007). In February 
2008, TeV flaring on ~ 1 day timescales was observed from M87, 
this time not associated with high activity from the knot HST- 1 . The 
correlation of TeV and X-ray fluxes on month to year timescales in 
the nucleus favors the nuclear region as the source of the TeV flares 
(Acciari et al. 2008; Harris et al. 2009). Even more convincingly, 
the 2008 TeV flaring was followed by new radio blobs moving out- 
ward from only ~ 100^ g from the black hole, suggesting that the 
gamma-ray flare comes from very close to the black hole (Accuari 
et al. 2009). 

If the jet of M87 is misaligned (pointing at an angle 9 ~ 30° 
away from the observer), the associated Doppler factor is 6 = 
1 /Tj ( 1 - VjCos9) ~ 1. In this case relativistic effects do not con- 
tribute to reduce variability timescales. The physical size of the TeV 
source is very compact, Z em < cf v ~ 5R g . This points to the possi- 
bility that the emission takes place in the vicinity of the black hole. 
TeV emission coming directly from the magnetosphere of the black 
hole is proposed by Neronov & Aharonian (2007; see also Levin- 
son 2000). A potential problem with this interpretation of the TeV 
emission is that infrared and optical radiation originating from the 
accretion disk may be strong enough to absorb TeV photons. The 
modeling of the emission of the accretion flow and fits to obser- 
vations of the nuclear emission indicate that TeV emission is ab- 
sorbed if produced at R < 20R g (Li et al. 2009). During the flaring, 
the high-energy emission extends up to at least ~ 10 TeV with the 
TeV spectrum well described by a power-law model with a photon- 
number index of T ~ -2.3 (Aharonian et al. 2006; Albert et al. 
2008). There is no evidence for absorption of the TeV emission, 
suggesting that the emission takes place at R cm > 2QR„. 

It has also been proposed that the TeV flares come directly 
from the inner (pc-scale) regions of a misaligned jet. These models 
are more complex than single-zone SSC models and either invoke 
jet deceleration on sub-pc scales (Georganopoulos et al. 2005) or 
a multi-zone configuration with a fast spine and slower outer layer 
(Tavecchio & Ghisellini 2008). One weakness of these models is 
that they tend to produce steeper TeV spectra than those observed, 
because of absorption by the dense soft (synchrotron) photon field 
along the line of sight. 

Our jets-in-a-jet picture provides a way out of these difficul- 
ties. While the jet is misaligned with respect to the observer, some 
minijets point outside the jet emission cone and in the direction to 
the observer. These "off-axis" minijets move relativistically toward 
the observer and can be compact, resulting in TeV flares with short 
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variability timescales. Furthermore, the relativistic beaming toward 
the observer allows the TeV emission to escape the production site 
without being significantly absorbed. 



3. 1 Jets in the jet of M87 

For more quantitative estimates we consider a jet with (isotropic) 
luminosity Lj iso that moves with bulk Tj. The jet is assumed to be 
strongly magnetized with a Poynting-to-kinetic flux ratio (mag- 
netization) cr » 1. As reference values, we use Tj = 5 and 
<x = 100. The isotropic equivalent jet luminosity is estimated to 
be Lj is0 = 4r j 2 Lj lrue ~ 10 45 erg/s. The energy density and magnetic 
field strength in the jet as function of radius are (e.g., Giannios et 
al. 2009): 

e\ = L i:im /4nr 2 cT 2 = 0.05/. ,<r. : l 5 2 erg/cm 3 , (7) 
and 

B\ = ^4tt7 = O.SI^r^irj Gauss, (8) 

where A = 10* A x and the spherical radius is R = r ■ R„ with i? g = 
4.5 x 10 14 cm, corresponding to the gravitational radius of a black 
hole of 3x 10 9 M o . 

We assume that a fraction of the magnetic energy of the jet is 
occasionally dissipated through reconnection. Our picture for rel- 
ativistic reconnection is the relativistic generalization of Petschek- 
type reconnection worked out by Lyubarsky (2005; see also Watan- 
abe & Yokoyama 2006; Zenitani, Hesse & Klimas 2009 for rela- 
tivistic MHD simulations that support this picture). High-cr mate- 
rial is advected into the reconnection region where the release of 
magnetic energy takes place. Part of the dissipated magnetic en- 
ergy serves to give bulk acceleration to a pair of minijets that move 
in opposite directions (in the rest frame of the jet) and the rest to 
accelerate particles in the minijets. We explore the possibility that 
emission from minijets results in the observed flaring activity. 

In this model, the material can leave the reconnection region 
with bulk r co close to the Alfven speed of the upstream plasma: 
r co ~ V^" — lOtr^' 2 in the rest frame of the jet (Petschek 1964; 
Blackman & Field 1994; Lyutikov & Uzdensky 2003; Lyubarsky 
2005). The energy density in the minijet is (Lyubarsky 2005) e em ~ 
ei = L 1Ji0 /4nr 2 cT 2 = 0.05/.. ^r. 'l ^ erg/cm 3 . 

Even though we consider a Poynting flux-dominated jet, the 
(downstream) emitting region is not necessarily magnetically dom- 
inated since a large part of the magnetic energy dissipates in the 
reconnection region. This has important implications for the re- 
sulting spectra. Actually, Lyubarsky (2005) shows that material 
leaves the reconnection region with r co ~ yTx for sufficiently 
small guide field. In this case, the magnetization of the minijet 
(downstream plasma) is cr cm = B 2 m /4ne em < 1, and the mag- 
netic field in the minijet rest frame is parameterized as 5 em = 
Vcr cm 4^e enl = O&crJ^L.'^r^Y^ Gauss. If the guide field of the 
upstream plasma is stronger, the minijet is strongly magnetized 
with cr em > 1, and slower, r co < \[tr. Furthermore, only a fraction 
~ 1/(1 + cr em ) of the magnetic energy is dissipated in the reconnec- 
tion region, leading to a rather weaker and slower minijets. Here 
we focus on fast minijets characterized by cr cm < 1. 



3.2 Particle distribution 

Assuming an electron-proton jet, the plasma coming into the recon- 
nection region (the upstream) contains ~ crm p c 2 magnetic energy 



per particle. Out of this energy a fraction yTx goes into bulk mo- 
tions with the rest being available to accelerate/heat particles. In 
the frame of the reconnection jet (minijet), the available energy per 
particle is ~ v7xm p c 2 . If electrons receive an appreciable fraction 
of the released magnetic energy / ~ 0.5, they are heated to charac- 
teristic thermal y factor 

Tech ~ / Vo-m p /m e ~ lO 4 /i/20-2 /2 , ( 9 ) 

assumed to be isotropic in the blob rest frame. 

In the context of the leptonic emission model discussed here, 
electrons of at least ~ 10 TeV (in the lab frame) are needed to po- 
tentially explain the observed emission that extends up to ~ 10 TeV. 
This means that the random component of the electron distribution 
should extend to y c > 2 x 10 7 /r em » y e ,ch ( m the minijet frame), 
where r em is the bulk Lorentz factor of the minijet. It is thus clear 
that electrons need to be accelerated well above the characteristic 
Lorentz factor y c ,ch- For simplicity, we will assume that the elec- 
trons follow a power-law distribution that extends to y c s> y ej ± . It 
remains to be shown that relativistic reconnection can lead to such 
a particle distribution. 

Since, for the parameters relevant to the jet of M87, particles 
with y c > y e ch are fast cooling (because of synchrotron and, mainly, 
inverse Compton emission; see next section), we consider emission 
coming from a rather steep particle distribution with p > 3. The 
particle distribution has a maximum cutoff y mllx . In our discussion, 
Tmax is not important (provided that it is high enough to produce 
~ 10 TeV photons through inverse Compton) and will be set to 
Tmax -» 00 ■ 

3.3 Radiation mechanisms 

The minijet contains energetic electrons that emit through syn- 
chrotron and synchrotron-self-Compton and, possibly, external in- 
verse Compton (EIC) mechanisms. We consider any source exter- 
nal to the minijet contributing seed photons that are scattered by 
the minijet as an EIC source. Following this definition, not only 
disk emission but also soft photons originating from other parts of 
the jet that interact with the minijet are labeled as EIC sources. 

In this Section, we explore the emission from a minijet for 
parameters of the model relevant to M87 (i.e., we set 6 = 30°, 
Tj = 5). Other parameters of the model (i.e., magnetization of the 
jet, radius of formation of the minijet) are less constrained. Since 
we model M87 as a misaligned blazar, for illustration we set these 
parameters to values inferred from modeling of Mrk 501 and PKS 
2155-304 (Giannios et al. 2009). The Lorentz factor of the minijet 
in the rest frame of the jet is, therefore, set to r co = yTr = 10. 
For these parameters, using eqs. (1), (2) we find that the minijet 
moves with r cm =12. The formation of the minijet is assumed to 
take place at R cm = 100/?„ and the magnetization of the plasma 
within the minijet is set to cr cm = 1/3 (though different values 
of cr and cr em are explored as well). From the observed luminos- 
ity (Lf ~ 10 42 erg/sec) and duration (7> ~ 10 5 sec) of the TeV 
flares, we can estimate the lab-frame energy contained in the mini- 
jet ~ LfTf/4T 2 m f ~ 5 x 10 44 erg and the typical length scale of 
the emitting region is / cm ~ {E cm IT aTt e a - n ) 117 ' ~ 10 15 cm. The mini- 
jet can be easily powered by the M87 jet of Lj >true ~ 10 43 erg/sec. In 
the following, we show that the minijet emits its energy efficiently 
in the ~ TeV energy range. 

The synchrotron emission of electrons with random Lorentz 
factor y ecil takes place at observed energy v syn ~ r cm y 2 v c = 
5Lj^a"2' 2 / 1 2 2 '"2 1 e V if the minijet points at us, where v c is the 
electron cyclotron frequency. The self-Compton emission appears 
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at v S sc ~ 7e v syn - O- 5 ^ '^/^^ ' GeV - Hi g her energy emis- 
sion is expected by particles with y c > y e ,min- 

In addition to the synchrotron emission in the minijet, ambient 
photon fields present in the nucleus of M87 are also a source of 
seed photons to be upscattered by the minijet. Some of these soft 
photons come from the accretion flow (Di Matteo et al. 2003; Li 
et al. 2009) and the rest from the inner jet. The observed nuclear 
emission has a flat Ef(E) spectrum in the range ~ 0.01 - 10 eV 
with integrated luminosity ~ 3 x 10 4 ' erg/sec declining at E < 0.01 
eV and E > 10 eV (e.g., Li et al. 2009). The nuclear emission has 
been resolved down to ~ 20if g in the radio (see, e.g., Spencer & 
Junor 1986; the radio emission is believed to come from the jet). 
The resolution is poorer in shorter wavelengths, corresponding to 
an upper limit to the size of the nucleus of ~ 100i? g in the mm and 
~ 3 x 10 4 tf g in the infrared (Why song & Antonucci 2004). 

The observed emission from the nucleus competes with the 
synchrotron emission from the minijet as sources of seed photons 
to be inverse-Compton scattered in the minijet. A conservative es- 
timate for the role of external soft photons comes from assuming 
that they are emitted isotropically from a region close to the black 
hole. If, on the other hand, a large fraction of the observed soft pho- 
tons comes from the jet (i.e. they are emitted isotropically in the jet 
frame), their intensity along the jet axis is much larger than that 
estimated assuming isotropic emission from the inner disk. 

If the soft photon field is emitted isotropically from a region 
with length scale R so i t , the typical angle 6 mX at which the soft pho- 
tons interact with the bulk motion of the minijet depends on the 
ratio Rioft/Rcm- For R so ft/R em ■< L the interaction angle equals ap- 
proximately the inclination of the minijet G mt ~ 6 ~ 30° while 
for Rsoft/Rem < 1, the interaction angle is larger, f int ~ 60°. The 
energy density of the external photon field at the location of the 
minijet is U ph ~ (1 - yS em cos e^dt^L^IAncR^ = 0.04(1 - 
j8 em cos #i„t)lt m l0 L S! 4 L5 r~ 2 erg/cm 3 (in the rest frame of the minijet), 
which appears to be of the same order of magnitude as that of the 
magnetic field, (?b = 8 x 10~ 3 Lj4 5 cr em ^r^Trl erg/cm 3 . It should 
be noted that the expressions for U p h and Ob hold for cr em < 1. For 
o"cm ~ 1 the magnetic energy density Ob saturates to a value that 
is independent of cr em , while for <x cm » 1, no minijet forms since 
the outflow from the reconnection region is expected to be slow 
and weak (since most of the magnetic energy is not dissipated by 
reconnection). 

For reasonable range of the parameters, either SSC or EIC can 
dominate the production of the gamma rays. First, we explore pa- 
rameters for which synchrotron seed photons and SSC dominate 
the bulk of the emission. We set cr em ~ 1 and assume a very com- 
pact emission region for the disk photons: -R so ft/^cm ■« 1- The bulk 
Lorentz factor of the minijet is taken to be rather modest, r cm = 8 
(by setting cr = 50). Having specified the properties of the minijet 
and of the ambient radiation field, we proceed with the calcula- 
tion of the resulting spectra. We calculate the synchrotron emission 
from a power-law distribution of relativistic electrons using stan- 
dard expressions. The SSC and EIC emission are calculated using 
the <5-function approximation for single electron emission. Klein- 
Nishina effects on the electron scattering are approximated by us- 
ing a step function for the energy dependence of the cross section: 
cr = cr r , for yjrvjm^c 1 < 3/4 and a = otherwise (see, for exam- 
ple Tavecchio, Maraschi & Ghisellini 1998 and Coppi & Blandford 
1990 for discussion on the accuracy of this method). 

The resulting spectra are shown in Fig. 4. Here we ignore any 
absorption of the TeV emission due to pair creation (discussed in 
the next section). The inverse Compton peak of the Ef(E) spectrum 
E p = r cn ,7 c _ c hm e c 2 corresponds to the inverse Compton emission 
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Figure 4. Synchrotron (dotted line) and inverse Compton (solid line) emis- 
sion from a minijet forming at R cm = 100i? e and beaming at an angle 
~ 30° with respect to the axis of the jet of M87. The magnetization of the 
minijet is set to tr cm = 1 and a power-law electron distribution is assumed 
with index p = 3.1. The observed emission from the nucleus is sketched 
with the dashed line and is assumed to come from a region R so [ t « R em . 
The black bow-tie reports the TeV spectrum observed during the early 2008 
flare (Albert et al. 2008). The EGRET 2cr upper limit is also shown (Sreeku- 
maretal. 1996). 



from electrons with the characteristic Lorentz factor y c-c h. Since the 
peak is located below the working range of the Cherenkov tele- 
scopes, E p < 100 GeV (Aharonian et al. 2006; Albert et al. 2008), 
we have the constraint that r em 7 c ch < 2 x 10 5 . Above the peak, the 
TeV spectrum is rather flat with photon-number index T = -2.3. 
The TeV spectrum depends on the index p of the electron distribu- 
tion (set to p = 3.1). The TeV emission becomes harder for lower 
p and vice versa. At E ~ 10 13 eV the spectrum steepens. Photons 
observed at E > 10 13 eV (i.e., E > 10' 2 eV in the rest frame of the 
minijet) mainly come from scattering in the Thomson regime of op- 
tical synchrotron photons (near-IR in the minijet frame), the energy 
density of which drops rapidly with decreasing energy. 

The slope of the TeV emission spectrum is compatible with 
that observed during the fast TeV flaring in 2005 (Aharonian et al. 
2006) and 2008 (Albert et al. 2008). At the same time the syn- 
chrotron emission is clearly pronounced in the optical, UV and 
X-rays, far exceeding (by a factor of 10-30) the typical nuclear 
emission. Note that, within the one-zone SSC model for the TeV 
emission discussed here, the synchrotron component peaks close to 
the optical and extends up to the hard X-rays rather independently 
of the adopted parameters. For photons to be upscattered to ener- 
gies E > 10 TeV in the Thomson regime, the synchrotron emission 
cannot peak much above the optical wavelengths. Furthermore, the 
electrons with random y c ~ 10 6 needed to produce the E > 10 TeV 
emission emit into the hard X-ray regime through synchrotron (for 
the magnetic field strength in the minijet of order 1 Gauss predicted 
by the model). Since powerful optical and X-ray flares with ~ 1 day 
timescales have not been observed in the nucleus despite the reg- 
ular monitoring of M87, we consider the SSC interpretation of the 
TeV emission unlikely. 

The synchrotron component of the minijet is weaker when 
cr cm < 1. This brings us to the parameter space for which the 
EIC mechanism plays a more important role. We set <r em = 1/3 
(or smaller) and assume an extended region for the external emis- 
sions: R so f t /R aTI < 1. The bulk Lorentz factor is set to r cm = 12 
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Figure 5. Same as in Fig. 4 but with the magnetization of the minijet set 
to o" em = 1/3 and the observed emission from the nucleus assumed to 
come from a larger region, i? so f t < R em . A power-law electron distribution 
is assumed with index p = 3.2. For these parameter the EIC mechanism 
dominates the gamma-ray emission. The green dash-dotted lines show the 
synchrotron and inverse Compton emission from an identical blob beam- 
ing within the jet axis as observed by an observer located at an angle 
9 = 1/Tj = 11° (with respect to the jet axis). 



(by setting cr = 100). The resulting spectra are shown in Fig. 5, 
where one can see that EIC dominates the GeV-TeV emission. For 
smaller values of cr cra the synchrotron and SSC components be- 
come weaker but the TeV emission remains practically the same. 
The TeV spectrum is rather well decribed with a power-law model 
with photon-number index T = -2.3. The slope of the TeV spec- 
trum depends on the spectrum of the soft seed photons (observed to 
be F so!t (E) oc ET 1 ) and the index p of the electron distribution (set 
to p = 3.2). At E ~ 3 x 10 14 eV the spectrum steepens because the 
energy density of the available soft seed photons with £ s < 10" 3 eV 
for scattering in the Thomson regime drops rapidly with decreasing 
energy. 

The slope of the TeV emission spectrum is compatible with 
that observed during the fast TeV flaring in 2005 (Aharonian et 
al. 2005) and 2008 (Albert et al. 2008). During the TeV flares the 
~ 100 MeV flux approaches the upper limits set by EGRET (see 
Fig. 5). If M87 is currently in a high GeV-TeV state, we predict 
that FERMI should soon have a significant detection of the nuclear 
emission. 

During flares, the synchrotron emission from the minijet dom- 
inates the nuclear emission in the X-rays, provided that <x em > 0.1. 
In this case, simultaneous X-ray flares of moderate strength are 
expected. Unfortunately, there have been no X-ray observations 
strictly simultaneous with the TeV flares. On the other hand, mon- 
itoring of the nucleus in X-rays in 2008 revealed that it was in a 
high X-ray state around the time of the TeV flares. Simultaneous 
multi-wavelength observations can greatly constrain magnetization 
of the emitting region. 

In Fig. 5, we also show for comparison the emission expected 
from M87 if it were observed at an inclination 6 = 1/Fj ~ 11°, 
i.e., at a = 1. In that case, the flares would be more powerful, 
powered by minijets emitting within the jet axis, and have observed 
durations of ~several hours. The synchrotron emission would peak 
in the extreme UV and the inverse Compton component at ~ 100 
GeV Both the timing and spectral properties would look similar to 
those of a low luminosity, high-frequency BL Lac (HBL) object. 



Note, however, that in HBLs the inverse Compton peak is typically 
not as pronounced with respect to the synchrotron one as that of 
Fig. 5. 



3. 3. 1 Escape of the Te V photons 

The spectrum of the TeV emission from M87 does not show any 
signature of absorption up to ~ 10 TeV. In this section, we examine 
the conditions under which the TeV emission from a misaligned 
minijet escapes pair creation in the region where it is produced. 

The dominant source of infrared-optical photons that can pair 
create with the TeV photons is the ambient nuclear emissiorf^. The 
size of the emitting region (which likely depends on the wave- 
length of the emission) and the location and direction of the minijet 
determine the typical angle (9j nt at which the TeV photons inter- 
act with soft photons. For a rough estimate of this angle we set 
^soft 5 femt resulting in mt ~ 60°. The pair-creation cross section 
peaks at a yy ~ cr T /5 for photons with total energy Eqm ~ 4m c c 2 
in the center of mass frame. The center-of-mass energy of a hard 
photon E h and a soft photon E s that form an angle (9 int ~ 60° is 
Eqm = V2£h£ s (l - cos 6) ~ -\JEhE s . The 77 annihilation cross 
section peaks at soft photon energy E s as 4/£'h,iTeV eV. Since the 
number density of soft photons increases at the lower energies, 
~10 TeV photons are the most likely to be absorbed by ~ 0.4 
eV target soft photons. The number density of soft photons is 
Marget = Lm/47rB% m cE s , where L iR ~ 10 41 erg/sec is the near- 
infrared luminosity of the nucleus (see Figs. 4, 5). Thus, the an- 
nihilation optical depth is 



l^IR,4l£h,10Tev/>2- 



(10) 



From this expression we conclude that for _R sofl < ^ cm , the ~ 10 
TeV photons marginally avoid significant attenuation if the minijet 
forms at i? e m ~ iOOR,,. If, however, R so f, is much smaller than R cm 
then escape is possible from smaller radii. 



4 DISCUSSION/CONCLUSIONS 

M87 is a well studied case of a misaligned AGN jet. Both the 
nucleus of M87 and the knot HST-1 have variable emission on 
timescales of weeks to months, in particular in the X-rays (e.g., 
Harris et al. 2009). More impressively, M87 reveals compact TeV 
emitting regions by producing flares on ~ 1 day timescales (Aharo- 
nian et al. 2006; Albert et al. 2008). VLBI observations during and 
after the 2008 TeV flares showed the ejection of two blobs within 
a distance of ~ 100^ g from the black hole, strongly suggesting the 
nucleus as the source of the flare (Acciari et al. 2009). 

The TeV flares may originate from the magnetosphere of the 
black hole (Levinson 2000; Neronov & Aharonian 2007) or a 
complex jet geometry (spine/layer interaction [Tavecchio & Ghis- 
ellini 2008]; collimating jet [Lenain et al. 2008]; decelerating jet 
[Georganopoulos et al. 2005]). Still, in models in which the jet is 
misaligned with respect to our line of sight, the Doppler factor of 
the emitting material is 6 ~ 1, implying the need to explain a very 
compact emitting region of length scale / cm < 5R g and the lack of 
TeV absorption. 



3 The synchrotron photons from the minijet may, in principle, also con- 
tribute to the TeV absorption. In practice, however, the strong beaming of 
the minijets toward the observer guarantees that the synchrotron photons 
can be ignored as a source of opacity. 
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Rapid gamma-ray flaring has been observed frequently in 
blazars. Mrk 501 and PKS 2155-304 show powerful flares at TeV 
energies with duration of several minutes, much shorter than the 
light crossing time of their central black holes (Aharonian 2007; 
Albert et al. 2007). This extremely fast flaring suggests that the 
flares do not reflect variability of the central engine in these jets 
but is connected to interactions within the jet, which lead to com- 
pact emitting regions and potentially relativistic motions in the rest 
frame of the jet. 

In this line of argument, we have previously proposed a jets- 
in-jet model for blazars (Giannios et al. 2009). Minijets, driven by 
magnetic reconnection and moving relativistically within the main 
jet, can power the fast evolving flares, make them transparent at 
TeV energies, and explain why the jet as a whole appears much 
slower. The reconnection may be triggered by MHD instabilities in 
the jet (like kinks) and/or reversals of the polarity of the magnetic 
field in the magnetosphere of the black hole or in the inner disk. 
When a field reversal takes place, interactions/collisions of parts of 
the jet with opposite polarity can lead to very efficient release of 
magnetic energy via magnetic reconnection. 

Reconnection within the jet produces two oppositely directed 
(in the jet's frame) minijets. One of them always points within the 
jet angle and is observable in blazars because their jets point at us. 
The other minijet (its counterpart) points outside the jet opening 
angle and is potentially observable to off-axis observers in case of 
misaligned jets, e.g., M87 or Cen A (Aharonian et al. 2009). 

The TeV flaring of M87 can be understood in this context. 
Furthermore, we predict that if M87 is currently in a high GeV- 
TeV state, FERMI should soon have a significant detection of the 
nuclear emission. Finally, if viewed on-axis, M87 would look like 
a fairly regular high peaked BL-Lac object with TeV (and possibly 
optical/Xray) flares of several hours. 

The high-energy emission from the minijets comes during 
their formation/ejection and/or the (short) cooling timescale asso- 
ciated with the TeV emitting particles. The observed TeV flares are 
likely part of a general increase of minijet events in the jet, most of 
which we do not see because they are beaming away from us. On 
a longer timescale, the minijets interact with the rest of the jet and 
are slowed down (with respect to the jet). During this interaction 
phase, shocks are expected to form, resulting in further accelera- 
tion of particles that can power an 'afterglow' emission. Interest- 
ingly the Acciari et al. (2009) VLBI observations that followed the 
February 2008 flaring of M87 showed the ejection of radio blobs 
from the nuclear region. It is tempting to associate these blobs with 
the collective afterglows from minijets brought to our view after 
their deceleration to the jet mean motion. 

By analogy with the ~ 5 min flares seen from Mrk 501 and 
PKS 2155-304, one might expect even faster evolving flares from 
M87 than the observed ~ 1 day ones. Taking into account that the 
black hole of M87 is a factor of ~ 3 larger than those of Mrk 501 
and PKS 2155-304 and that the large inclination of the M87 jet 
dilutes timescales by a factor of a 1 ~ 10, one could still expect 
flares evolving in as short as 2 - 3 hours from M87. 
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